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■^1- i Abstract. A spinning particle in the Schwarzschild spacetime deviates from 

geodesic behavior because of its spin. A spinlcss particle also deviates from 
geodesic behavior when a test radiation field is superimposed on the Schwarzschild 
background: in fact the interaction with the radiation field, i.e., the absorption 
O , and re-emission of radiation, leads to a friction-like drag force responsible for 

^J , the well known effect which exists already in Newtonian gravity, the Poynting- 

Robertson effect. Here the Poynting-Robertson effect is extended to the case of 
spinning particles by modifying the Mathisson-Papapetrou model describing the 
motion of spinning test particles to account for the contribution of the radiation 
force. The resulting equations are numerically integrated and some typical orbits 
are shown in comparison with the spinless case. Furthermore, the interplay 
^. ■ between spin and radiation forces is discussed by analyzing the deviation from 

circular geodesic motion on the equatorial plane when the contribution due to the 
radiation can also be treated as a small perturbation. Finally the estimate of the 
amount of radial variation from the geodesic radius is shown to be measurable in 
principle. 
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1. Introduction 

The motion of classical spinning test particles in a given gravitational background 
is described by the well known Mathisson-Papapetrou (MP) model [1] [2]. Let 
U a — dx a /dr be the timelike unit tangent vector to the "center of mass line" Cjj 
of the spinning particle used to perform the multipole reduction, parametrized by the 
proper time r. The equations of motion are 

-T- = - ^u a pU v S al3 = F^ in ^ , (1.1) 

dr Z 

— — = P^U" - P U U» , (1.2) 

dr 

where P M is the total 4-momentum of the particle and the antisymmetric tensor S^ v 

denotes the spin tensor (intrinsic angular momentum) associated with it; both fields 

are defined only along this center of mass world line. This system of 10 equations 

evolves P and S along Cjj but contains 13 unknown quantities: U (3), P (4), S (6). 
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Consistency of the model is ensured by imposing the Tulczyjew-Dixon [3J HI HJ [8] 
supplementary conditions 

ST„ = 0. (1.3) 

Moreover, implicit in the model is the requirement that the spin structure of the 
particle should produce very small deviations from geodesic motion, in the sense that 
the length scale naturally associated with the spin should be very small when compared 
with the one associated with the curvature tensor of the spacetime itself. Otherwise, 
large values of spin would require taking into account the particle backreaction on the 
spacetime metric, i.e., the problem should be approached from a completely different 
point of view. 

Let us consider a spinless test particle orbiting a star which emits radiation. The 
radiation pressure of the light emitted by the star, in addition to the direct effect of 
the outward radial force, exerts a drag force on the particle's motion. This usually 
causes the body to fall into the star unless it is so small that the radiation pressure 
pushes it away from the star, a phenomenon called the Poynting- Robertson effect since 
it was first investigated by Poynting 9 using Newtonian gravity and then calculated 
in the framework of linearized general relativity by Robertson |10| . Successively many 
authors studied the Poynting-Robertson effect in more concrete situations, starting 
from the case of slowly evolving elliptical orbits for meteors [IT] , to more recent works 
[T2l [T3l HH1 [T51 1TB] . where rotation of the emitting star is taken into account. More 
recently the Poynting-Robertson effect for a spinless test particle orbiting a black hole 
was studied in both the Schwarzschild and Kerr spacetimes 171 without the restriction 
of slow motion, but ignoring the finite size of the radiating body. 

Here we generalize the above discussion to the more realistic case of a spinning 
test particle subject to the Poynting-Robertson effect by including the radiation forces 
in the Mathisson-Papapetrou model. Today we have evidence of the existence of 
accreting matter around massive compact objects, e.g., active galactic nuclei [18] . 
The dynamical behavior of particles in close orbits around massive objects, while 
interacting with their radiation field, can be relevant when studying the evolution of 
shell or disk-like configurations of dust around intense radiative relativistic sources, 
where the loss of angular momentum via the Poynting-Robertson effect could act 
as a dust accretion mechanism. The radiation mechanism around a real accreting 
compact object is generally very complicated. We will limit ourselves to the case 
of a coherent flux of photons traveling along geodesies in some preferred direction. 
Possible scenarios include a hot neutron star, a black hole accreting radiation or a 
system with an accretion disk which radiates. When the approximation of point-like 
test particles is no longer valid, we expect that the real Poynting-Robertson effect 
offers a different behavior with respect to the standard Poynting-Robertson effect for 
small dust particles. The Mathisson-Papapetrou model allows us to take into account 
the actual size of the particle in the framework of general relativity by introducing a 
characteristic length of the particle itself through its spin. 

It is worth mentioning that there also exists a wide literature concerning pseudo- 
classical test spinning particles, whose equations of motion reduce under certain 
limit to the classical MP equations [Tj|J [13 US HI]- In fact, spinning particles can 
be equivalently described by pseudo-classical mechanics models in which the spin 
degrees of freedom are characterized in terms of anticommuting Grassmann variables, 
associated — in the semiclassical limit — with the components of the spin tensor of the 
particle. The Lagrangian formulation can be used as well to study spinning particle 
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motion in external fields. Here, however, we consider only classical test bodies, leaving 
for future work further generalizations of this analysis to the case of pseudo-classical 
particles. 

2. Motion in the Schwarzschild spacetime 

Consider a Schwarzschild spacetime, whose line element written in standard 
coordinates is given by 

ds 2 = -iV 2 di 2 + N~ 2 dr 2 + r 2 (d9 2 + sin 2 dd(/) 2 ) , (2.1) 

where N = (1 — 2M/r) 1 ' 2 denotes the lapse function, and introduce the usual 
orthonormal frame adapted to the static observers (or Zero Angular Momentum 
Observers, ZAMOs) following the time lines 

&i = N~ l d t , e f = Nd r , e* = -d s , ei = — ^A , (2.2) 

° r v rsva.0 

where {<9t, d r , 8$, d^,} is the coordinate frame. 

We limit our analysis to the equatorial plane 9 = n/2, where the situation 

is relatively receptive to analytical treatment. As a convention, the physical 

(orthonormal) component along — ds which is perpendicular to the equatorial plane 

will be referred to as "along the positive z-axis" and will be indicated by the index I, 

when convenient: ei = — e&. 

2.1. Test particles subject to PR effect 

Let a pure electromagnetic radiation field be superposed as a test field on the 
gravitational background described by the metric (|2.1[) . with the energy-momentum 
tensor 

T af} = <5> 2 k a k p , k a k a =0, (2.3) 

where k is assumed to be tangent to an affinely parametrized outgoing null geodesic 
in the equatorial plane, i.e., k a \/ a k^ = with k e = 0. 

We will only consider photons in the equatorial plane which are in outward radial 
motion with respect to the ZAMOs, namely with 4-momentum 

k = E(ri)(n + e?) , (2.4) 

where n = e t ~ is the ZAMO 4- velocity and E(n) — E/N is the relative energy of the 
photon as seen by the ZAMOs. Here E = —k t is the conserved energy associated with 
the timelike Killing vector held dt- Note also that L = k^ = 0, i.e., the conserved 
angular momentum associated with the rotational Killing vector field 9^ is assumed 
to vanish. This is consistent with having a non-rotating light source in a non-rotating 
spacetime. 

Since k is completely determined, the coordinate dependence of the quantity $ 
then follows from the conservation equations T a ^ -^ — 0, and will only depend on r in 
the equatorial plane due to the axial symmetry. We find 

$ = ^. (2.5) 

r 

Consider now a test particle moving in the equatorial plane 6 = 7r/2 accelerated 

by the radiation field, i.e. with 4-velocity 

U = j{U, n) [n + v(U, n)\ , v{ U, n) = v r e? + i/*e; = ^(sin ae? + cos ae i ) , (2.6) 
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where 7(17, n) = 1/yl — \W(U, n)\\ 2 = 7 is the Lorentz factor and the abbreviated 
notation v a = v(U,n) a has been used. Similarly v = \\v(U,n)\\ and a are the 
magnitude of the spatial velocity v(U, n) and its polar angle measured clockwise from 
the positive <p direction in the r-(f> tangent plane respectively, while v = v{U, n) is the 
associated unit vector. Note that a — corresponds to azimuthal motion with respect 
to the ZAMOs, while a = ±7r/2 corresponds to (outward/inward) radial motion with 
respect to the ZAMOs. 

A straightforward calculation gives the coordinate components of U 

* dt 7 dr * a dd> 71^ 

where r is the proper time parameter along Cjj, and U e = d9/dr — 0. Solving these 
for the magnitude and polar angle leads to 



tana = _L^!: v= iJ(*)\m(¥)\ (2.8) 

Nr d<f) ' N 2 ]] \dtj \dtj y J 

The scattering of radiation as well as the momentum-transfer cross section a 
(assumed to be a constant) of the particle are assumed to be independent of the 
direction and frequency of the radiation so that the associated force is given by 

nana nsi 

F ( - rad '> a = -aV(U) a pT^U^, (2.9) 

where V(U) a p = <52 + U a Up projects orthogonally to U. Explicitly 

F (rad)a = - a &(p(JJ)'* p tf) (k^) , (2.10) 

implying 

^ (l ' ad) = ^2 7 3 (1 - "*) [(?* ~ "> + (l-v f - (>) 2 ) e? - (1 - u f )Je } ] , (2.11) 

where we have used the notation ctQqE 2 — mA, as in Ref. [17]. Test particle motion 
is then described by the equation 

mailjy = m— — = F (rad) ^ , (2.12) 

dr 

and has been studied in detail in Ref. |17j . 

2.2. Generalization to spinning particles 

The most direct and simple generalization of Eq. (|2.12[) to the case of spinning test 
particles consists in including the radiation force term in Eq. ()l.lj) . so that one has 

T)PA« 

= ^(spin)^ _|_ ^(rad)p (2.13) 

dr 
plus the additional relations (|1 . 2|) and (|1.3[) involving the spin. Let us proceed to 
analyse the motion of spinning particles subject to the Poynting-Robertson effect in 
the equatorial plane of Schwarzschild spacetime. 

The 4-momentum P — mu for motion in the equatorial plane is 

u = 7u[n + p u (sma u er + cosa u ei)} , j u = — ^== , (2.14) 

V 1 ~ K 
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and introduce the spin vector associated with S^ by spatial duality 

S? = Ucff^Sp, , (2.15) 

where r] a p 7 s = ^/— ge a fi 1 s is the unit volume 4- form and e a p 7 s (eoi23 = 1) is the 
Levi-Civita alternating symbol. It is also useful to consider the scalar invariant 

s 2 = \S» V S» V , (2.16) 



constant along Cu because of Eqs. (jl.21) and fll.3[) . Consistency of the model requires 
that the length scale |s|/m associated with the spinning particle be much smaller than 
the one associated with the background spacetime, say M , namely 

'" «1. (2.17) 



1 ' mM 

Let us consider Eqs. (|1.1[) and (|1.2[) with U a given by Eq. (|2.7|) and u a given by 
Eq. (|2.14[) . In the spinless case P is aligned with U, i.e., u = U, implying that v = v u . 
The presence of the spin causes a change in both U and u according to 

U = Uo+sU s , u = U + su s , (2.18) 

where 

U =*f (n + vfa + v$e}) (2.19) 

satisfies Eq. (J2.12J) and corrections are first order in the spin. Higher order terms in 
Eqs. (jl.ip and (|1 .2[) are neglected. This leads to two different sets of equations for 
zeroth and first order in spin respectively, which are listed in | Appendix A[ 
The spin force to first order in s is given by 

F(spin) = _3mM^. 7 2^ ( ^ n + &f) (2 2Q) 

r o 
We find that the mass of the spinning particle to is a constant of motion. Furthermore, 
from the evolution equations for the spin it follows that the spin vector has a single 
nonvanishing and constant component along 9 (or z), namely 

S = -S S e s = se £ . (2.21) 

Figures [TH3] show some numerical solutions for the orbits in the strong field region. 
Of course it only makes sense to consider the exterior solutions for radii larger than 
some minimum radius R outside the horizon in order to model the geometry outside 
a star (or some other physical source) of radius R producing the outflow of radiation. 

In the case of spinless particles there exists a condition representing the balancing 
of the gravitational attraction and the radiation pressure at constant ro and 4>q, namely 

A / 2M\ 1/2 1M 



M-\}-—) ^ r " = r ^ " 1 - AVM' ■ (2 - 22) 

This behavior also characterizes the motion of spinning particles as well, as shown in 
|Appendix A.l| Figs. Q] and [2] show some typical solution curves starting initially with 
purely azimuthal velocity either inside (Fig. [IJ or outside (Fig. [5]) the critical radius 
at which a particle initially at rest with respect to the ZAMOs (which in turn are 
at rest with respect to the coordinate system) remains at rest. For comparison, the 
corresponding curves for spinless particles with the same initial data are also shown. 
If A/M -C 1, the critical radius approaches the horizon r( cr it) « 2M. For instance, 
an initially Keplerian circular orbit gradually spirals towards the central source, as 
illustrated in Fig. [3l 
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Figure 1. The orbit of a spinning particle (solid curve) subject to the Poynting- 
Robertson effect is shown for the choice of parameters A/M = 0.8 and s = 0.5 
(X = r cos <f> and Y = rsin0 are Cartesian-like coordinates). The starting point 
is located at r (0) = 4Af and <j> (0) = with u u0 (0) = 0.7, a u0 (0) = 0, t s (0) = 0, 

r s (0) = and cf> s (0) = 0, i>J(0) = and uf(0) = 0. The values of the spin 
parameter has been exaggerated in order to distinguish the difference from the 
motion of a spinless particle (dashed curve). The inner circle is the horizon 
r = 2M, while the outer circle is at the critical radius n crit ) = 5.5M which 
is outside the initial data position. 



3. Deviation from the circular geodesic 

Consider now the corrections to geodesic circular motion, by taking the effect of the 
radiation field to also be small. 

In the absence of both spin and radiation we assume the geodesic motion of the 
particle to be circular at r = r$ (tq > 3M in order Uk to be timelike), that is 

U = U K = 7ic(n±^e^), (3.1) 

where the Keplerian value of speed (vk) and the associated Lorentz factor (7^) and 
angular velocity (Of) are given by 



AI 



r - 2M 



V K = \ TTT, 1K = \ TTT' Ck = \~3- (3-2) 

V ro — 2M V ro — SM y rg 

The ± signs in Eq. (|3.ip correspond to co-rotating (+) or counter-rotating (— ) orbits 
with respect to increasing values of the azimuthal coordinate <$> (counter-clockwise 
motion as seen from above). The azimuthal direction in the local rest space of Uk 
pointing in the direction of relative motion (i.e., the boost of e 1 in the local rest space 
of Uk) is specified by the following unit vector orthogonal to Uk in the t-<p plane 

U K =lK{v K n±e^), (3.3) 

where the ± signs are correlated with those in Uk- 
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Figure 2. The orbit of a spinning particle (solid curve) subject to the Poynting- 
Robertson effect is shown for the choice of parameters A/M = 0.6 and s = 0.5. 
The starting point is located at ro(0) = AM and </>o(0) = with v u q(0) = 0.5, 

a u0 (0) = 0, t s (0) = 0, r s (0) = and S (O) = 0, i/j(0) = and uf(0) = 0. 
The corresponding orbit for a spinless particle is also shown (dashed curve). The 



critical radius r 



(crit) 



3.125M is inside the initial data position. 
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Figure 3. The orbit of a spinning particle (solid curve) subject to the Poynting- 
Robertson effect is shown for the choice of parameters A/M = 0.01 and s = 0.5. 
The starting point is located at ro(0) = AM and (f>o(0) = with ^ u o(0) = vk ~ 

0.7071, q„o(0) = 0, t s (0) = 0, r s (0) = and <t> s (0) = 0, <(0) = and vf(0) = 0. 
The corresponding orbit for a spinless particle is also shown (dashed curve). In 



this case r 



(crit) 



2M. 
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The parametric equations of Uk are 

tic = to + —t = to + Tkt , 

a 7r 
1KVK 



b K = O ± -^^t = O ± n K r , (3.4) 

T"0 



where now to, tq and </>o are constants and 



r K = .H F °~, n K = ±.J - . (3.5) 

V r o - 3M ' r V r - 3M V ; 

It is convenient to introduce a friction parameter /, so that the length scale A 
associated with the radiation held is much smaller than M, i.e., 

Therefore, in the present analysis corrections to geodesic motion will be limited to 
first order terms in both parameters s and /, according to 

t = t K + ftf + sts , r = r + fr/ + sr§ , <f> = <f>K + f4>f + s4>s , 

v* = /«/> + |«/> , J = ±i/ x + /i/* + szi , 

i/„ = ± v K + fu uf + §is uS , a u = fa uf + sa uS , (3.7) 

where tx and 4>k are given by Eq. (|3.4p . This implies 

U = U K + fU f + SUs , (3.8) 

where 

Uf = f -i/k — ± 7^/ ) U K + lKv)er , 

U s = (-VK V f ± T^-ff ) t/if + 7^£e f . (3.9) 

Similarly 

u = U + fu f + su s , (3.10) 

with [/ given by Eq. (|3. 81) and 

as discussed in | Appendix B| 

To first order in s and / the spin force and radiation force are given by 

F< s P in > = T 2mM$r? K <? K v K et , 

F( rad ) = - mfil K v K ('YKVKUK ~ et ) , (3.12) 

respectively. The ratio between the magnitudes of these forces has the form 

] C^ = 4(-) V \f^ E - (3 . 13 ) 

Its behavior as a function of ro in units of \s\/ f is shown in Fig. |4] 
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Figure 4. The ratio between the magnitudes of spin and radiation forces given 
by Eq. 113.131 1 is plotted in units of \s\/f as a function of rg/M. 



The equations governing first order perturbations are listed in |Appendix B| The 
corresponding solution is given by 

MP K 



"• = T -» n. 



• sin(Oe P r) , 



'0 "op 



'K 



O 



K 



sin(O cp T) + 2r a ( K — [cos(O ep T) -I}} , 



■3 

K is K [cos(n cp T) - 1], 



3MCJ 

— — i/K[<MS\,itepi 



"op 
,3 






XI 



, , [cos(0 C pT) - 1] - ^ sin(O cp r) + Vl\t 



(3.14) 



and 



6M 2 ^3 
is = T p-^- sin(S2 cp T) - S2 cp T , 



"ep 
13 



t, = 4r C^i§f |[cos(O cp r) - 1] + ^-^-[ 
"o P I <Zr <, K ll K 



sin(n ep r) - n cp r] + -~f 2 K n 2 cp T 2 } , 



rg = ±3r — 2 — [cos(S2 cp r) - 1J , 

^o P 

r f = - r o(K~jTf- < [cos(O cp t) - 1] - 2r ( K — ^[sin(O cp r) - O cp r] \ , 
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± ^-t s , 0/ = ±%i/ , (3.15) 



V K V K 



where 



IM(r -6M) 
rg(r -3A/) 



^ = 1-371 777T ( 3 - 16 ) 



is the well known epicyclic frequency governing the radial perturbations of circular 
geodesies. 

The constant term in r$ represents the slight change in the radius of the circular 
orbit about which the solution oscillates with proper period 27r/51 cp . In contrast, the 
presence of a secular term in r/ is responsible for the deviation from geodesic motion 
due to friction, which is measurable in principle. In fact, by taking the mean values 
over a period of the perturbed radius we can estimate the amount of variation of the 
radial distance 

\ r I r UL 



1 - aTrroC*^ ] / ^/-/^kAu 

* 'cp 



"ep 

For instance, for the motion of the Earth about the Sun we find 



(3.17) 



(*L\ « / T 2 x io-" (f/!!!k « 3 x 10- 15 T 4 x 10" 15 « 10" 15 , (3.18) 

\ r I cm 

since r « 1.5 x 10 13 cm, M — M Q w 1.5 x 10 5 cm and the ratio (s/m)® ~ 200 cm for 

the Earth; the friction parameter is related to the ratio between the solar luminosity 

£ « 3.8 x 10 33 erg/s and the Eddington luminosity QUEEZ] £ Ed d ~ L3 x l ° 38 er §/ s ' 

and for the Sun is given by / w 3 x 10 -15 . Therefore, in this case the effect of the 

radiation field on the orbit is of the same order as that due to spin. Note that the 

estimate of the contribution due to spin is in agreement with [23] . 

The effect of the spin may become important when the orbiting extended body 

is a fast rotating object. To illustrate the order of magnitude of the effect, we 

may consider the binary pulsar system PSR J0737-3039 as orbiting Sgr A*, the 

supermassive (M ~ 10 6 M@) black hole located at the Galactic Center [241E5], at a 

distance of r ~ 10 9 Km. The PSR J0737-3039 system consists of two close neutron 

stars (their separation is only dAB ~ 8 x 10 5 Km) of comparable masses tua — 1.4 M Q , 

raj ~ 1.2 Mq), but very different intrinsic spin period (23 ms of pulsar A vs 2.8 s of 

pulsar B) 26 . Its orbital period is about 2.4 hours, the smallest yet known for such an 

object. Since the intrinsic rotations are negligible with respect to the orbital period, 

we can treat the binary system as a single object with reduced mass \iab — 0.7 M & 

and intrinsic rotation equal to the orbital period. The spin parameter thus turns out 

to be equal to s w 1.0 x 10~ 3 . The luminosity of Sgr A* is about 10 3 £q, whereas its 

Eddington luminosity is £Edd ~ 10 n £©, so that / m 10 -18 . Therefore, in this case 

' 5r\ 
r I 

Therefore, in this case the effect of the spin on the orbit dominates with respect to 
the friction due to the radiation field. 

4. Concluding remarks 

We have studied the motion of a classical spinning body in the field of a central 
radiating object. The model adopted is the standard Mathisson-Papapetrou model 



. =a 7.6 x 10" 19 T 1.8 x 10~ 7 « 10~ 7 . (3.19) 
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suitably modified by accounting for the contribution of the Poynting-Robertson 
radiation force in the equations of motion. We have numerically integrated the 
whole set of Mathisson-Papapetrou equations in the case of equatorial motion and 
coherent flux composed of radially emitted photons. The spin vector turns out to 
have only a constant nonvanishing component orthogonal to the motion plane. We 
have shown some typical solution orbits in comparison with the spinless case. The 
latter is characterized by the existence of a critical radius at which the balancing of 
the gravitational attraction and the radiation pressure occurs at constant radial and 
azimuthal coordinates depending on the strenght of the radiation field. This feature 
has been proved to be maintained also in presence of spin. Dust particles would 
congregate at this radius leading to rings of matter to form. 

Furthermore, we have discussed the interplay between spin and radiation forces by 
analyzing the deviation from circular geodesic motion on the equatorial plane when 
also the contribution due to friction can be treated as a small perturbation. The 
features of the motion thus depend on two different parameters, the spin parameter 
and the friction parameter, which are taken as small in order to avoid backreaction. 
The presence of the spin causes a slight change in the radius of the circular orbit 
about which the solution oscillates and an increase/decrease of the angular velocity 
depending on whether the particle is co/counter rotating (i.e., moving clockwise or 
anticlockwise with respect to the positive <j) direction, respectively). In contrast, the 
presence of a secular term in the radial deviation due to friction determines a spiraling 
behavior of the orbit. This leads to a (average) radial variation from the geodesic 
radius whose amount is measurable, at least in principle. 

The model presented here allows to account for the finite size of the particle 
subject to the Poynting-Robertson effect in a framework which is genuinely relativistic. 
Obviously, in order to be physically more realistic, the model should be further 
generalized to take into account, for instance, the finite size of the radiating source and 
the contribution of higher order multipoles in the description of the actual size of the 
orbiting body and its shape, e.g., by including quadrupolar terms in the Mathisson- 
Papapetrou equations of motion. 

Appendix A. Solving the MP equations: general case 

The parametric equations for the center of mass line C\j are given by 

t = t + st s , r = r + sr s , <f> = </> + s<j) S , (A.l) 

so that 

u f = v r +sut, v $ = ui + §u{, (A.2) 

where all quantities are functions of the proper time r. A similar expansion holds for 
u, i.e., 

v u = v u0 + sv uS , a u = a u0 + sa u § , (A.3) 

where is u o — v§. The first order correction to u turns out to be 

luoVuoVusn + {-v r Q a U s + cos a U QV uS )e? + (tft a uS + sin a u0 v uS )e7 , 



7«o 



InoViUo) 



(-"Ja«s + cosa„ ^s)ef + {vq a uS + s'ma u0 iy U s)e l j > , (A.4) 



where V(Uq) projects orthogonally to Uq 
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The zeroth order quantities satisfy Eq. (J2.12I) . i.e., 
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dr 

dv U Q 

dr 
da M o 

dr 



No 



C.KVR 



cos a u0 

VuQ 



dr 
dr 
sin a u0 



7«o^o^o 



»o _ 7«o^o 



dr 



'•() 



TT — (1 - ^)( sina «0 - I'uo) 

M ro 



r/K M r 



where 



(A.5) 



v'o = Vv& sin a u0 , 1$ = u u0 cos a u0 , (A. 6) 

and the Keplerian value of speed vk and the associated Lorentz factor jk an d angular 
velocity (k have been introduced in Eq. (|3.2I) . 

The hrst order quantities satisfy the equations 
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(A.7) 
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The remaining quantities u u g and a u $ are related to the first order spatial velocities 
by the algebraic relations 

AM cosa„o 



ToiVo 7«o 
AM 7 u0 
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rgiVo U u0 
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V U Q 



(A.5 



Appendix A.l. Equilibrium solutions 

In order to find an equilibrium position at a given point (t"(t#), tt/2, 4>{t*)) for values 
t > t* of the proper time we have to impose first the conditions dr/dr = and 

d0/dr = 0, which are fulfilled by Vq = = z^f (i.e., i/ u0 = 0) and v r ~ = = fj . 
Requiring their first derivatives with respect to r to be identically vanishing as well 
gives the condition (J2.22I) from the zeroth order equations (|A.5|) . The first order set 
(|A.7|t yields 

,2 
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(A.9) 




The last condition is nothing but Eq. (|2.22[) . whereas the previous one implies r 
at all values of the proper time. Equilibrium positions are thus characterized by 

U = n , P = m{n + MsC, K v K e~^) , (A.10) 



at r = r 



(crit) ■ 



Appendix B. Solving the MP equations: the case of small / 

The quantities first order in s satisfy Eqs. (|A.7[) and (IA.8I) where also terms 
proportional to fs are neglected, taking into account that 



Vq = , vl = ±v K 
according to Eq. p.7[) . We have 
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(B.2) 
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The remaining quantities u u g and a u $ are related to the first order spatial velocities 
by the algebraic relations 



vk 



(B.3) 



The equations for the quantities first order in / come from the linearization of 
Eq. (IA.5I) . where also terms proportional to fs are neglected 
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Finally, the remaining quantities v u f and a u / turn out to be simply given by 



(B.4) 



'-V 



a«/ = ±- 



7 



(B.5) 



We have now to solve the coupled system of equations (|B.2[) and (|B.4|) . Taking 
the derivative of the equations for the first order radial components of the spatial 
velocity with respect to proper time yields 



d',J 



2 ,,f 



,2 /-2 



dr 2 



K v vl = o , -^ + fi^i/j; + 2r^^a = o 



(B, 



where the epicyclic frequency f2( ep ) has been introduced in Eq. (|3.16|l . The general 
solution of Eq. (|B.6|) is straightforward 



ai sin(il cp r) + 02 cos(51 cp r) , 



7 



61 sin(f2 cp r) + b 2 cos(f2 cp r) - 2T K v K 



2 ^K 



(B.7) 



cp 



The first order azimuthal components of the spatial velocity are then given by 
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Finally, the first order corrections to the orbit turn out to be 
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* "cp 



-[a 2 sin(r2 cp r) — a\ cos(Q, cv ,t)] + c\ , 



rf = — [6 2 sin(riop'r) - b\ cos(^ cp t)] - 2 — —j-t + c 2 . 

* 'Cp ^*0 ^ 'oTl 



"op 

± ^T^i sin (^c P r) + a 2 cos(r! op T)] T T^w^i T MN u k ]t + e x , 



2 -z4f[ai sin(fi cp r) + a 2 cos(tt cp r)] T |^| 

2 7|Cx r , . . , n S , . , n X, 7l l; ,,, , ■ v . 



</>/ = ± -^§p[6i sin(fi cp r) + 6 2 cos(fi cp r)] T ^73 [Af + 3 7if ^^c 2 ]r 

r 0"ep ^ r iV 



±^oCl§fr 2 + e2 , (B.9) 

where 01, a 2 , 61, 6 2 , ci, c 2 , dj., d 2 , ei, e 2 are arbitrary integration constants. Their 
values are fixed by requiring that all first order quantities vanish at r = 0. The 
corresponding solution is given by Eqs. (I3.14J) and (|3.15[) . 
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